Astronomy & Astrophysics manuscript no. 60Fe'wwang 


February 1, 2008 


(DOI: will be inserted by hand later) 





SPI observations of the diffuse 60 Fe emission in the Galaxy 

S 

W. Wang 1 , M.J. Harris 1 , R. Diehl 1 , H. Halloin 2 , B. Cordier 4 , A.W. Strong 1 , K. Kretschmer 1 , J. 
. Knodlseder 3 , P. Jean 3 , G.G. Lichti 1 , J. P. Roques 3 , S. Schanne 4 , A. von Kienlin 1 , G. Weidenspointner 3 , C. 

Wunderer 5 

, 1 Max-Planck-Institut fur extraterrestrische Physik, Postfach 1603, 85740 Garching, Germany 

2 AstroParticule et Cosmologie (APC), 11 Place Marcelin Berthelot, 75231 Paris, France 

3 Centre d'Etude Spatiale des Rayonnements, B.P. 4346, 31028 Toulouse Cedex 4, France 
DSM/DAPNIA/SAp, CEA Saclay, 91191 Gif-sur-Yvette, France 
Space Sciences Laboratory, UC Berkeley, CA 94720, USA 



CO ■ 4 



43 

0- Received 

o 

' Abstract. Aims. Gamma-ray line emission from radioactive decay of Fe provides constraints on nucleosynthesis 

' in massive stars and supernovae. 

Methods. The spectrometer SPI on board INTEGRAL has accumulated nearly three years of data on gamma-ray 
emission from the Galactic plane. We have analyzed these data with suitable instrumental-background models 
and sky distributions to produce high-resolution spectra of Galactic emission. 

Results. We detect the 7-ray lines from 60 Fe decay at 1173 and 1333 keV, obtaining an improvement over our 
\jr^ ' earlier measurement of both lines with now 4.9 a significance for the combination of the two lines. The average 

^\ \ flux per line is (4.4 ± 0.9) x 10 -5 ph cm -2 s -1 rad -1 for the inner Galaxy region. Deriving the Galactic 26 Al 

gamma-ray line flux with using the same set of observations and analysis method, we determine the flux ratio of 
CO ; 60 Fe/ 26 Al gamma-rays as 0.148 ± 0.06. 

, Conclusions The current theoretical predictions are consistent with our result. We discuss the implications of 

these results for the widely-held hypothesis that 60 Fe is synthesized in core-collapse supernovae, and also for the 
' closely-related question of the precise origin of 26 Al in massive stars. 

o ' 
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X. 

^ , 1. Introduction et al. 1992). With its much longer half-life (~ 10 6 years), 

■ - - 26 Al may propagate over significant distances of ~ few 

The radioactive isotope 60 Fe is believed to be synthesized hundred pC; &nd accumulates in the interstellar medium 

through successive neutron captures on Fe isotopes (e.g. (IgM) from many supernovae; until i njec tion and /3-decay 

5b Fe) in a neutron-rich environment inside massive stars, &re ^ b&lance ^ the IgM) giying rige tQ ft diffuse and eyen 

before or during their final evolution to core collapse su- Galaxy _ wide glow (Mahoney et al . 1982; p ra ntzos and 

pernovae (CCSN). Due to its long decay time (r ^2.2 My), DieU igg6) _ The behavior q{ 60 p e from the same masgive _ 

faU Fe survives to be detected after the supernova ejected it ^ sQmces shoM foUow ^ q{ 26a1j gince itg half life ig 

into the interstellar medium, by /3-decay via &0 Co and 7 similarj ^ 2 .2 x 10 6 years. In addition, 60 Fe could also be 

emission at ,1173 keV and 1333 1 keV - like other radioactive produced in substantial amounts by rare subtypes of SN 

isotopes: 44 Ti, 56 ' 57 Co, and 26 Al. These isotopes provide I& (Woosley 1997)j which would then be point sources of 

evidence that nucleosynthesis is ongoing in the Galaxy 60p e g amma _ ra y S 
(with CCSN occurring at intervals ~ 30 — 100 yr, The et 

al. 2006; Diehl et al. 2006). The discoveries of the above Although the detections of these isotopes are in agree- 
isotopes have fulfilled predictions of nucleosynthesis the- ment with the broad outlines of CCSN nucleosynthesis 
ory, at least in its general outline. The 44 Ti and 56 ' 57 Co theory, there are discrepancies in the details. For exam- 
isotopes have relatively short half-lives (< 100 years) and pie, 44 Ti lives long enough so that it should have been 
are thus detected as point-like sources, e.g. in young Type detected from several recent Galactic supernovae if these 
II supernova remnants (SNRs) Cas A and 1987A, respec- occur at a rate of ~ 2 per century; but with the excep- 
tively (Iyudin et al. 1994, 1995; Matz et al. 1988; Kurfess tion of Cas A (Iyudin et al. 1994, Renaud et al. 2004, 

Renaud et al. 2006), no other SNRs have been detected 
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al. 1999, Pliischke et al. 2001) confirms that it follows 
the overall Galactic massive star population, but concen- 
trations of 26 Al are seen in very young OB associations 
where even the most massive stars are not expected to 
have exploded (e.g. Cygnus OB associations). This has 
led to the suggestion that most 26 Al production is asso- 
ciated with ejections in an earlier, pre-explosion phase of 
stellar evolution (Knodlseder et al. 1999). 

Measurements of 60 Fe promise to provide new infor- 
mation about these issues, specifically the massive star 
nucleosynthesis in the late pre-supernova stages. RHESSI 
reported observations of the gamma-ray lines from 60 Fe 
with 2.6 g significance, and an average flux of (3.6 ± 
1.4) x 10~ 5 ph cm~ 2 s- 1 (Smith 2004). The analysis of the 
first year of data from the SPI spectrometer on ESA's 
INTEGRAL spacecraft resulted in a similarly marginally- 
significant detection of these 7-ray lines from 60 Fe (~ 3<r, 
Harris et al. 2005), with an average line flux of (3.7 ± 
1.1) x 10~ 5 ph cm" 2 s" 1 . SPI on INTEGRAL has accu- 
mulated more data since then; here we analyzed 2.5 years 
of data, aiming at a consolidation of the INTEGRAL/SPl 
measurement of 60 Fe gamma-rays. 

In the following section, we describe our observations 
with SPI. In §3, we describe the steps and methods of 
data analysis. The scientific results, including the spectra 
and intensities of both 1173 keV and 1333 keV lines, and 
the flux ratio of 60 Fe/ 26 Al are presented in §4. Discussions 
and a summary follow in §5. 

2. Observations 

The INTEGRAL spacecraft was launched on October 17, 
2002, into a high-inclination, high-eccentricity orbit in- 
tended to avoid the increased background from the Earth's 
trapped radiation belts. INTEGRAL'S orbital period is 3 
days. The spectrometer SPI consists of 19 Ge detectors 
actively shielded by a BGO anti-coincidence shield. It has 
a tungsten coded mask in its aperture which allows imag- 
ing at ~ 3° resolution within a 16° x 16° full coded field 
of view (imaging on INTEGRAL is mainly performed at 
lower energies by the IBIS telescope, with which SPI is co- 
aligned). The Ge detectors are sensitive to gamma-rays 
between 15 keV and 8 MeV, with a total effective area 
~ 70 cm 2 at 1 MeV, and achieve an energy resolution of 
~ 2.5 keV at 1 MeV (Roques et al. 2003, Attie et al. 2003). 
However, cosmic-ray (CR) impacts degrade this resolution 
over time, and the instrument is periodically shut off for 
a few days while annealing (by heating from cryogenic 
temperatures to 100°C) is applied to the detectors to re- 
store the energy resolution by thermal curing heating of 
the CR-induced defects (Roques et al. 2003, Leleux et al. 
2003). 

In space operations, INTEGRAL with its IBIS and 
SPI telescopes is pointed at predesignated targets, with 
a fixed orientation for intervals of typically ~ 2000 s (re- 
ferred to as pointings), which are successively arranged as 
a standard pattern of neighbouring pointings ~ 2° apart 



(dithering), and covering target region of interest for im- 
proved imaging (Courvoisier et al. 2003). 

The interval covered by the observations analyzed here 
is December 3, 2002 — August 22, 2005 (orbits 17-359). 
The instrument was in operation for most of this time, 
but important gaps were caused by five of the annealing 
episodes described above, and by the regular perigee pas- 
sages with gaps due to the radiation belts. The live time 
was further reduced by the failure of two of the 19 detec- 
tors during this interval (December 2003 and July 2004) 
and by the occurrences of solar flares. 

3. Data analysis 

The basic measurement of SPI consists of event messages 
per photon triggering the Ge detector camera. We distin- 
guish events which trigger a single Ge detector element 
only (hereafter single event, SE), and events which trigger 
two Ge detector elements nearly simultaneously (hereafter 
multiple event, ME), which is the case for ~ 30% of 60 Fe 
gamma-ray photons. 

The analysis steps are: (1) Assembly of data selected 
to be free of contamination by, e.g., solar-flare events. (2) 
Modelling the instrumental background. (3) Fitting the 
measured data in narrow energy bins with the background 
model and a model of celestial gamma-ray emission, folded 
through the instrumental response into the data space 
of the measurement; the fitted amplitude of the celestial 
model per energy bin then comprises the spectrum of ob- 
served sky emission. (4) Deriving line parameters of the 
celestial signal. From the two event types, and from the 
two gamma-ray lines of 60 Fe, we obtain four independent 
measurements of 60 Fe emission. 

In the following, we discuss the analysis steps and 
their implementation in the corresponding software utili- 
ties within MPE's INTEGRAL data analysis system. 

3.1. Data selection and assembly 

We test data quality per each 'science window' (typically a 
time interval of ~ 30 minutes corresponding to one point- 
ing), applying selection limits to 'science housekeeping' 
parameters such as the count rates in several onboard ra- 
diation detectors, instrument status codes, data owner- 
ship, and orbit phase (utility spi select sew). We employ 
the INTEGRAL Radiation Environment Monitor (IREM; 
Hajdas et al. 2003), the SPI plastic scintillator anticoinci- 
dence counter (PS AC), and the rate of saturating events 
in SPI's Ge detectors (from events depositing > 8 MeV 
in a detector; hereafter referred to as GEDSAT rates) to 
exclude solar-flare events and other erratic background 
increases. Regular perigee background increases are addi- 
tionally eliminated through a 0.05-0.95 window on orbital 
phase. From the selected events, spectra are accumulated 
per each detector and pointing, and together with dead 
time and pointing information assembled into the analy- 
sis database. 
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We establish databases of spectra per pointing per de- 
tector for single and multiple events (SE and ME), in 
spectral ranges with 1 keV bin size in ~ 20 keV bands 
around each of the 60 Fe lines, and in adjacent energy bands 
for determination of instrumental background (see below). 
Typical such raw data spectra around two 7-ray lines are 
presented in Figure 1, and show that instrumental back- 
ground including strong instrumental lines dominates the 
raw spectra. The energy ranges around the two lines ex- 
tend from 1153 - 1193 keV and 1313 - 1359 keV. For the 
background modelling from the adjacent continuum, we 
use 1163 - 1169 keV plus 1177 - 1184 keV, and 1318 - 
1328 keV plus 1336 - 1349 keV energy bands, respectively. 
From orbits 17 to 359, a total observation time of ~ 24 
Ms is thus obtained, with a data set consisting of 14623 
pointings with one spectrum per each detector and event 
type. 

3.2. Background modelling 

SPI spectra are dominated by the intense background ra- 
diation characteristic of space platforms undergoing CR 
bombardment. Much of this radiation is prompt, result- 
ing directly from CR impacts, whose variation with time 
ought to follow that of the incident CR flux. Other compo- 
nents arise from radioactive isotopes produced by the CR 
impacts, whose decay lifetimes are long compared to the 
coincidence window of the detector-anticoincidence elec- 
tronics. Radioactive isotopes will increase in abundance 
until CR production balances /3-decay. A time series of 
their 7-ray emission will be the convolution of the time 
behavior of the prompt CR source with an exponential 
decrease from decay. Local radioactivity in the spacecraft 
and instruments themselves thus will generate both broad 
continuum background emission and narrow gamma-ray 
lines from long-lived radio-isotopes. Varying with energy, 
background components may exhibit complex time vari- 
ability due to their origins from more than one physical 
source. 

We derive models for the background contribution per 
energy bin and pointing in each detector from compar- 
isons to presumed 'tracers' of background; for these, we 
use independent contemporaneous measurements aboard 
INTEGRAL (see Halloin et al. 2007, Diehl et al. 2006, 
also background model studies by Sizun et al. 2007). 
Comparing different candidate tracers to observed back- 
ground variations, we find the best tracer(s) of background 
for each energy band. Candidate tracers with high sta- 
tistical precision are the rates measured in SPI's plas- 
tic scintillator anticoincidence detector, the rate of sat- 
urated events in the BGO shield, the detector-by-detector 
GEDSAT rates (from events depositing >8 MeV in a Ge 
detector), and Ge-detector rates integrated over rather 
wide energy bands (for better statistical precision). For 
our analysis, we make use of the GEDSAT rates tracing 
the prompt CR activation to the 60 Fe lines (see above). 



Each such tracer is regarded as a 'template' for time 
variability of background. By normalizing such a template 
to the set of counts per detector spectrum (i.e. per point- 
ing), a background model can be constructed, which is ap- 
plicable to the analysis data set of spectra. We construct 
our background model from two components, one derived 
from adjacent energy bands, and an additional component 
modelling instrumental line contributions. 

The most important radioactive background line com- 
ponent for our analysis is the emission from 60 Co decay in- 
side the instrument, whose two lines at 1173 and 1333 keV 
are actually part of the celestial 60 Fe decay chain we want 
to measure. The convolution of the GEDSAT CR source 
tracer with an exponential function e~*' T , in which the 
decay time r is that of 60 Co (7.6 yr), is expected to give 
a good background model for energy bins containing con- 
tributions from these lines, i.e. around 1173 and 1333 keV 
(see Figure 2). Another possible radioactive contaminant 
in our spectral region of interest may be the strong 69 Ge 
K-shell electron capture line (1337 keV, Wcidcnspointner 
et al. 2003) which blends into the 1333 keV 60 Fe line; its 
lifetime is 2.35 days. Because of the short decay time is 
short, we just take the GEDSAT as a background compo- 
nent tracer for this strong instrumental line feature. 

Thus, in a first step, the detector-by-detector count 
rates in the four continuum bands are fitted by the 
GEDSAT time series to construct an 'adjacent-energies' 
background template. Then, the set of detector-by- 
detector spectra per pointing in 1 keV bins covering the 
~ 20 keV intervals around and including the 60 Fe lines are 
fitted to the sum of adjacent continuum template, plus 
the radioactivity template for 60 Co, plus an additional 
GEDSAT template to capture any additional prompt 
background. Thus we generate a background model for 
our actual set of spectra, which is constrained in its vari- 
ability by the various background tracers. Contributions 
of celestial gamma-ray events are small, their variability is 
intrinsically different due to the coded-mask modulation 
of the background model supplemented by dithering, and 
thus is expected only to affect global normalization, but 
not the variability with time. 

We then orthogonalize the different background model 
components for an improved convergence of the fitting 
(also see §2.2 of Diehl et al. 2006; this step ensures that 
each higher-order background component carries new and 
independent information ). 

Representing our data sets with background-models 
only, we obtain reduced \ 2 values of 1.185 and 1.194 for 
the 1173 and 1332 keV line bands (SE, 277824 d.o.f), and 
0.665 and 0.663 for the 1173 and 1333 keV line bands (ME, 
615038 d.o.f.), respectively. Note that for the low num- 
ber of counts in ME spectra, x 2 statistics does not apply. 
Therefore, in our model fitting approach (§3.3), we use 
Poissonion statistics and maximize the Poisson-likclihood 
function. 
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Fig. 1. Raw data spectra around the energies of two 60 Fe lines in one-INTEGRAL-orbit observations (3 days), rep- 
resenting the instrumental lines and continuum background. For the 1173 keV case (left), three strong instrumental 
lines are obvious: 44 Sc (1157 keV), 60 Co (1172.9 keV), 182 Ta (189.4 keV), and the 60 Co line blends with the 60 Fe line 
(Weidenspointner et al. 2003). For the 1333 keV case (right), the 60 Co background line (1332.5 keV) and the other 
strong instrumental line of 69 Ge (1336.8 keV) blend with the 60 Fe line, an instrumental line at 1347 keV also comes 
from 69 Ge electron captures. 
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Fig. 2. Background coefficients for the 60 Co radioactivity build-up model around the 1173 keV (left) and 1333 keV 
(right) lines, respectively. 




Fig. 3. Two plausible sky distribution models adopted for diffuse 60 Fe emission from the Galaxy. Both have been 
derived from COMPTEL 26 Al measurements (Pliischke et al., 2001; Knodlseder et al., 1999). The image obtained 
with the Maximum-entropy method (left) shows more structure, while the image derived with the Multi-resolution 
Expectation Maximization Method (MREM) appears smooth from its noise-filtering, accepting only image structure 
which is enforced by the measurements. 



3.3. Model fitting and spectra 

We combine above background models with a spatial 
model for sky emission (Figure 3) to fit our data, allowing 
for adjustments of fit parameters for background and sky 
intensities. For this, we use a maximum-likelihood method 
(spimodfit, for more details of methods, see Strong et al. 
2005) based on Poissonion statistics. We this derive per 
energy bin, the fitted parameter values with uncertain- 
ties, the covariance matrices, and the fitted model com- 



ponents. The counts per energy bin, per detector, and per 
pointing are fitted to the background model described in 
§3.2 and the assumed sky map of celestial emission (e.g., 
26 Al distribution obtained by COMPTEL, see Figure 3) 
as convolved into the domain of spectra per detector, en- 
ergy bin, and pointing through the pointing sequence and 
the instrumental response: 

D e , d , p = £ A^JJ m ' n + £ ^ B U P + Se, d , P , (1) 

m,n i—l 
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Fig. 4. Residuals of counts after fitting, in different projections. The upper figure shows the residuals versus time 
(days of the Julian Date, starts at 1 Jan 2000), where counts per pointing have been re-binned into 3-day intervals for 
clarity (here for single events SE in the 1163 — 1182 keV energy band). The lower two figures present the residuals 
projected with the energy (SE cases for the bands around the 1173 and 1333 keV lines separately). Two different time 
periods are shown (dashed: first half of database; solid: second half of database). For the 1333 keV line case, because of 
the strong instrumental 69 Ge line at 1337 keV, there still exist the large residuals around 1337 keV. Residuals around 
zero generally confirm that our background models are adequate except for the instrumental line around 1337 keV, 
with x 2 /d-o.f. in model fittings (§3.3): 1.182 (SE) and 0.670 (ME) for the 1173 keV line; 1.191 (SE) and 0.667 (ME) 
for the 1333 keV line. 



where e,d,p are indices for data space dimensions: en- 
ergy, detector, pointing; m,n indices for the sky dimensions 
(galactic longitude, latitude); A is the instrument response 
matrix, / is the intensity per pixel on the sky. Coefficients 
S for the sky map intensity (constant in time) and 
Pb,2, 0b,3 for three background intensities are derived (see 
§3.2 and Figure 2; different normalizations allowed for 
each camera configuration of 19/18/17 functional detec- 
tor elements due to failure of two Ge detectors) , and the 
amplitude /3 S comprises the resultant spectra of the sig- 
nal from the sky. S is the count residue after the fitting. 
Generally, a good fit will lead to residuals being statisti- 
cally distributed around zero (see Figure 4). 

Separate fits are made for each of the Fe lines and for 
SE and ME, respectively. These produce four spectra of 



sky brightness amplitudes per energy bin. If background 
models are adequate, the continuum energy bands should 
show sky amplitudes statistically consistent with zero sky 
brightness. Non-zero sky brightness should show up in the 
form of lines with a shape conforming to the instrumen- 
tal line width, possibly broader if the celestial line energy 
is broadened already (as had been discussed for 26 Al, see 
Diehl et al. 2006). Due to the low intensity of expected 
60 Fe emission, we do not attempt to derive line shape in- 
formation at all, and rather derive the 60 Fe line ampli- 
tude by fitting a Gaussian with fixed instrumental width 
to the spimodfit result spectra, determining line centroids 
and amplitudes in this way. From strong nearby instru- 
mental lines, we determine that a Gaussian with a width 
of 2.76 keV is an adequate model for narrow lines in the 
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regime of the 60 Fe lines; we use this shape to determine 
60 Fe line parameters from the spimodfit result spectra. For 
a consistent celestial signal from 60 Fe, the lines should be 
around the 60 Fe decay energies of 1173 and 1333 keV and 
equal (within errors) in amplitudes for all four data sets. 

4. Results 

4.1. SPI detection of m Fe 

From independently fitting the SE and ME data sets at 
1 keV-wide energy bins in the energy ranges of the two 
lines of 60 Fe emission (1173 keV and 1333 keV), we obtain 
the four spectra from spimodfit (§3.3) shown in Figure 5. 
Our fitted model consists of the sky intensity distribution 
of 26 Al from 9-year COMPTEL observations (Figure 3, 
maximum entropy image as a standard input, Pluschke et 
al. 2001), together with the background model based on 
measurements in energy bands adjacent to the 60 Fe lines 
and on background tracers (see above). We find excess 
amplitudes for the celestial component in each of the four 
spectra. 

Residual background imperfections are small ex- 
cept for the instrumental line at 1333 keV. The 
continuum flux levels are close to zero within 
10 _5 ph cm~ 2 s _1 rad -1 keV -1 , without apparent 
energy dependence, and adding a linear slope to our 
spectral fit function does not improve the fit . The 69 Ge 
line at 1337 keV has also been eliminated rather well, 
though not completely Suppression is a factor ~ 4 better 
than our previous work (Harris et al. 2005). However, this 
strong instrumental line still leads to a spectral signal 
which is recognized clearly as an artifact and thus not 
confused with the celestial line. 

We then determine the 60 Fe line parameters in these 
spectra fitting a continuum and Gaussians representing 
the SPI spectral response (§2). The width of the Gaussian 
profiles has been fixed at 2.76 keV, the value derived 
from simultaneously-measured adjacent instrumental- line 
shapes (e.g. the instrumental line at 1107 keV). We thus 
obtain 60 Fe line positions and intensities from fitting the 
four spectra of Figure 5, where the total line flux for each 
line in the inner Galaxy region (e.g. —30° < I < 30° and 
— 10° < b < 10°) is determined from the intensity of our 
fitted Gaussians and from the normalization of the input 
sky map. We present line fluxes of each of the two 60 Fe 
lines for the SE and ME data in Table 1, respectively. All 
line flux values are consistent within uncertainties. 

A superposition of the four spectra of Fig. 5 is shown in 
Fig. 6. Line energies of the 60 Fe lines in the laboratory are 
1173.23 and 1332.49 keV. For this superposition, we there- 
fore define the zero of the relative energy axis at 1173 and 

The diffuse gamma-ray continuum emission around 
60 Fe line energies in the inner Galaxy is ~ 2 x 
10" 6 ph cm -2 s _1 rad -1 keV" 1 from COMPTEL measure- 
ments (Strong et al. 1999), which is well below the error bars 
of spectral counts in Figure 5. 



1333 keV, to derive the summed spectrum of all 60 Fe sig- 
nals. The line flux estimated from the combined spectrum 
is (4.4 ± 0.9) x 10 _5 ph cm -2 s^ 1 rad -1 . Our significance 
estimate for the combined spectrum is ~ 4.9<r, adding un- 
certainties of the individual spectra in quadrature. This 
improves upon earlier 60 Fe signal reports from RHESSI 
(Smith 2004) and the first year of SPI data (Harris et al. 
2005). 

The signal of celestial 60 Fe is very weak, and 
marginally significant (< 3<r) in each of our four spectra. 
Therefore we cannot evaluate line shape information. The 
lines in our spectra appear well represented by Gaussians 
with instrumental widths (Fig. 5 and 6), suggesting that 
the celestial 60 Fe lines are intrinsically-narrow lines. If we 
assume a line broadening of 1 keV, line fluxes would in- 
crease by ~ 6%. 

4.2. 60 Fe and 26 Al 

It has been argued that 26 Al and 60 Fe share at least some 
of the same production sites, i.e. massive stars and super- 
novae (e.g., Timmes et al. 1995; Limongi & Chieffl 2006). 
In addition both are long-lived radioactive isotopes, so 
we have good reasons to believe their gamma-ray distri- 
butions are similar as well. Therefore we adopt the sky 
distribution of 26 Al gamma-rays as our model for celestial 
60 Fe gamma-ray distribution. 

From our model-fitting approach (§3.3), in principle, 
the spectral result depends on the input sky map (also 
see Eq. 1). We do not know the real distribution of 60 Fe 
sources in the Galaxy, but we need a realistic model for 
these sources to derive correct spectra and line intensity. 
As an alternative to the COMPTEL 26 Al maximum en- 
tropy and MREM sky maps, we also tried a different 
sky map: an exponential-disk model with scale radius 4 
kpc, and scale height 180 pc. This should also be a good 
first-order representation of 26 Al emission from the inner 
Galaxy, and but avoids the fine structure inputs from the 
COMPTEL image measurements which may partly arise 
from instrumental or analysis imperfections. We find that 
different input sky maps do not change the line profiles 
and intensities significantly (variations of the line fluxes 
are within 10% of the quoted values, i.e. below the un- 
certainties). If we use very different and scientifically im- 
plausible sky distribution models, e.g., a point source at 
I = ±20° , a bulge model (a Gaussian profile in the 
Galactic center), or a COMPTEL 26 Al map zeroed in the 
inner region of -40° < I < 40°, -10° < b < 10°, then no 
signals are found near the 60 Fe line energies. 

If we hypothesize a massive-star origin for 60Fe, then 
the measurement of the gamma-ray flux ratio 60 Fe/ 26 Al is 
important for discussions of the astrophysical nucleosyn- 
thesis origins of the two radioactive isotopes, and the nu- 
clear physics involved in models for their production (ad- 
dressing the uncertain nuclear reaction cross sections and 
half-lives). For this purpose, we apply the same analysis 
method on the 26 Al data from the same observations, an- 
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Fig. 5. The spectra of two gamma-ray lines of 60 Fe from the inner Galaxy: 1173 keV and 1333 keV. We have shown 
the results both from SE and ME databases. The data points are fitted with Gaussian profiles of fixed instrumental 
width (2.76 keV), and fixed continuum slope (flat), with x 2 / rf -°-/- of nts: 1-06 (SE) and 1.09 (ME) for the 1173 
keV line and 1.16 (SE) and 1.11 (ME) for the 1333 keV line. For the SE database, we find a line flux of (4.2 ± 
1.6) x lCT 5 ph cm" 2 s" 1 rad" 1 for the 1173 keV line and (3.5 ± 1.5) x l(T 5 ph cm" 2 s" 1 rad" 1 for the 1333 keV line. 
For the ME database, the line flux is (5.8 ± 1.9) x lCT 5 ph cm" 2 s" 1 rad" 1 for the 1173 keV line and (5.2 ± 2.1) x 
l(T 5 ph cm" 2 s' 1 rad" 1 for the 1333 keV line (also see Table 1). 



alyzing the 1785 — 1826 keV energy band (1 keV bins). 
We again generate a specific background model (from ad- 
jacent energy bands and the GEDSAT background inten- 
sity tracer), and apply the same input sky map (the 9- 
year COMPTEL 26 Al maximum entropy map) in model 
fitting the SE data. This yields a flux ratio of 60 Fe/ 
26 Al in a self-consistent way, with F( 60 Fe)/F( 26 Al) = 
(14.8 ± 6.0)%. Here, the uncertainty has been estimated 
from the respective model-fitting uncertainties of the two 
SE databases. Alternatively, from the combined spectrum 
of 60 Fe lines (Fig. 6), and adopting the 26 Al intensity 
measured with SPI before (,F( 26 A1) = (3.04 ± 0.31) x 
10~ 4 ph cm~ 2 s _1 rad -1 for the inner Galaxy, see Fig. 
5 of Diehl ct al. 2006), we obtain F( 60 Fe)/F( 26 Al) = 
(14.5 ±4.0)%. 

4.3. The Cygnus and Vela regions 

The Cygnus region is one of the most active nearby star 
formation regions in our Galaxy. It contains a large num- 



ber of massive stars and rich OB associations at a distance 
of 1-2 kpc (e.g. Knodlseder 2000; Le Duigou & Knodlseder 
2002), and has been recognized as a prominent source of 
26 Al with COMPTEL (Diehl et al. 1995a, Pluschke et al. 
2001) and with INTEGRAL/SPl (Knodlseder et al. 2004). 
From these massive stars, 60 Fe gamma-rays at the 1173 
keV and 1333 keV decay lines may also be expected. Since 
the majority of Cygnus region star-clusters arc young (~ 3 
Myr), from population synthesis studies of the massive 
stars in the Cygnus region it has been suggested that the 
60 Fe production is low, consistent with the small num- 
ber of recent supernova events inferred for this region 
(F(1173 keV) - 2 x 10" 6 ph cm" 2 s" 1 , see Knodlseder 
et al. 2002). 

The Vela region in the southern sky includes even 
more nearby massive stars, and prominent candidate 26 Al 
sources such as the Wolf-Rayet star WR 11 in the bi- 
nary system 7 2 Vel at d ~ 260 pc (van der Hucht et al. 
1997), and the core-collapse supernova remnant Vela SNR 
at d ~ 250 pc (Cha et al. 1999). 1809 keV line cmis- 
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Fig. 6. The combined spectrum of the Fe signal in the inner Galaxy, superimposing the four spectra of Figure 2. In 
the laboratory, the line energies are 1173.23 and 1332.49 keV; here superimposed bins are zero at 1173 and 1333 keV. 
We find a detection significance of 4.9o\ The solid line represents a fitted Gaussian profile of fixed instrumental width 
(2.76 keV), and a flat continuum. The average line flux is estimated as (4.4 ± 0.9) x 10 _5 ph cm" 2 s" 1 rad -1 . 




Fig. 7. Flux ratio of the gamma-ray lines from the two long-lived radioactive isotopes 60 Fe/ 26 Al from several observa- 
tions, including our SPI result (also see Table 2), with upper limits shown at 2a for all reported values, and comparison 
with the recent theoretical estimates (the upper hatched region from Prantzos 2004; the horizontal line taken from 
Timmes et al. 1995; the lower hatched region, see Limongi & Chieffi 2006). Our present work finds the line flux ratio 
to be (14.8 ± 6.0)%. Note that the primary instrument results on 60 Fe is compromised in the figure by different and 
uncertain 26 Al fluxes used for the 60 Fe/ 26 Al ratio (see Table 2 for details). 



sion of 26 Al from the Vela region has been reported from 
COMPTEL (Diehl et al. 1995b). And this emission is not 
attributed to the nearby 7 2 Vel source (Oberlack et al. 
2000) and may be either diffuse or due to supernovae in 
the region. Therefore the Vela region should be a good 
candidate for gamma-ray line emission from 60 Fe. 

Splitting our model sky map into independent com- 
ponents to fit the data, we do not see significant contri- 



butions from either region; our estimated upper limits for 
60 Fe gamma-rays from the Cygnus (68° < I < 92°, -10° < 
b < 14°) and Vela (255° < I < 273°, -9° < b < 9°) re- 
gions are ~ 1.1 x 10~ 5 ph cm" 2 s" 1 (2a). More detailed 
studies are in progress (Schanne et al. 2007, and Martin 
et al., in preparation). 
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5. Summary and discussions 

We report the detection of 60 Fe decay gamma-ray lines 
in the Galaxy from 2.5 years of SPI observations. With 
our new measurements we detect both the 1173 keV 
and 1333 keV line of 60 Fe from SPI single and multiple- 
detector events. Combining our four spectra from inde- 
pendent model fits we obtain an 60 Fe signal from the 
Galaxy with a significance of 4.9a. This improves upon 
our previous measurements (see Harris et al. 2005). The 
average 60 Fe line flux from the inner Galaxy region is 
(4.4 ± 0.9) x 10~ 5 ph cm~ 2 s _1 rad -1 , assuming intrinsi- 
cally narrow lines and a sky distribution equal to that 
of 26 Al as measured by COMPTEL (Phischke et al. 
2001). From the same observations and analysis proce- 
dure applied to 26 Al, we find a flux ratio of 60 Fc/ 26 Al 
of (14.8 ± 6.0)%. The 60 Fe signals are too weak to deter- 
mine line shape details; it appears that a Gaussian with 
the width of the instrumental resolution can fit the data 
well. This would imply that broadening of 60 Fe lines from 
astrophysical processes is not significant; most of the 60 Fe 
may be distributed in a rather normal interstellar medium 
with turbulent velocities below ^300 km s _1 . In order to 
investigate the variations over the Galaxy, we search for 
60 Fe emission from the Cygnus and Vela regions, and do 
not find 60 Fe signals. 

Many experiments and efforts were made (see Table 
2 and Fig. 7) to measure the 60 Fe gamma-ray emis- 
sion that was predicted by theory - we now provide the 
most significant detection to date. In Fig. 7, we show the 
previous constraints on the flux ratio of 60 Fe/ 26 Al to- 
gether with the result of this work, and compare the ob- 
servational results with different theoretical predictions. 
The earliest observational limit was given from HEAO- 
3, F( 60 Fe)/F( 26 Al) = 0.09 ± 0.08, an upper limit be- 
ing 0.27 (Mahoney et al. 1982) (in Fig. 7, we chose to 
give limits at 2a for all reported values below a signifi- 
cance of 3cr). Another limit was obtained with the SMM 
Gamma- Ray Spectrometer, a flux ratio of 0.1 ± 0.08, the 
upper limit being - 0.27 (Leising & Share 1994). OSSE 
aboard the COMPTON Observatory gave a flux ratio 
of 0.21 ± 0.15, and the upper limit is ~ 0.51 (Harris 
et al. 1994, 1997). COMPTEL aboard the COMPTON 
Observatory also found 60 Fe gamma-rays, and reported a 
flux ratio value of 60 Fc/ 26 Al of 0.17 ± 0.135, which trans- 
lates into an upper limit ~ 0.44 (Diehl et al. 1997). The 
Gamma-Ray Imaging Spectrometer (GRIS) reported an 
upper limit for the ratio of < 0.14 (2a, Naya et al. 1998). 
RHESSI then reported the first detection of 60 Fe gamma- 
ray lines, and gave a flux ratio 0.17 ± 0.13 with two-year 
data (Smith 2004). The first year data of SPI gave a flux 
ratio 0.11 ±0.07 (Harris et al. 2005), and the present anal- 
ysis of the 2.5-year SPI data finds a flux ratio 0.148 ±0.06. 

Theoretical predictions of the ratio of 60 Fe/ 26 Al have 
undergone some changes since Timmes et al. (1995) pub- 
lished the first detailed theoretical prediction. In their pa- 
per, they combine a model for 26 Al and 60 Fe nucleosyn- 
thesis in supernova explosions with a model of chemical 



evolution, to predict that the steady production rates are 
(2.0±1.0)M Q Myr- 1 for 26 Al, and (O.75±O.4)M Myr" 1 
for Fe, which corresponds to a gamma-ray flux ratio 
F( 60 Fe)/F( 26 Al) = 0.16 ± 0.12. This prediction would 
be consistent with our present measurements. Since 2002, 
theoreticians have improved various aspects of the stellar- 
evolution models, including improved stellar wind models 
and the corresponding mass loss effects on stellar struc- 
ture and evolution, of mixing effects from rotation, and 
also updated nuclear cross sections in the nucleosynthe- 
sis parts of the models. As a result, predicted flux ratios 
60 Fe/ 26 Al rather fell into the range 0.8±0.4 (see Prantzos 
2004, based on, e.g. Rauscher et al. 2002, Limongi & 
Chieffi 2003) - such high values would be inconsistent 
with several observational limits and our SPI result (but 
see Woosley & Heger 2007 for comments on nuclear re- 
action rate updates ). Recently, new calculations of 26 Al 
and 60 Fe yields from massive stars of solar metallicity have 
been presented (Limongi & Chieffi 2006; Heger & Woosley 
2007). Limongi & Chieffi (2006) combined their individ- 
ual yields, using a standard stellar-mass distribution func- 
tion, to produce an estimate of the 60 Fe/ 26 Al gamma-ray 
flux ratio expected from massive stars. Their calculations 
yield again a lower prediction for the 60 Fe/ 26 Al flux ra- 
tio of 0.185 ± 0.0625, which is again consistent with the 
observational constraints (see Figure 7). 

In summary, uncertainties still exist, both in models 
and measurements of 60 Fe. On the theory side, stellar evo- 
lution in late stages is complex, nuclear reactions include 
neutron capture on unstable Fe isotopes, and explosive 
nucleosynthesis adds yet another complex ingredient. On 
the experimental side, cosmic ray induced 60 Co radioac- 
tivity in the instrument and spacecraft and the limita- 
tions of spatial resolutions and sensitivity are issues re- 
flected in the substantial uncertainties in experimental 
values. With more INTEGRAL/SPl data to come, and 
also with the development of next-generation gamma-ray 
spectrometers/telescopes, gamma-ray observations hope- 
fully can help with an independent view on the astrophys- 
ical model components. 
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